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We study experimentally a system comprised of linear chains of spin- 1/2 nuclei that provides a 
test-bed for multi-body dynamics and quantum information processing. This system is a paradigm 
for a new class of quantum information devices that can perform particular tasks even without 
universal control of the whole quantum system. We investigate the extent of control achievable on 
the system with current experimental apparatus and methods to gain information on the system 
state, when full tomography is not possible and in any case highly inefficient. 



I. INTRODUCTION 

The dynamics and control of complex quantum many- 
body systems has elicited renewed interest within the 
field of quantum information science. Not only is there a 
need to improve coherent control of an increasing number 
of qubits to build a scalable quantum computer, but also 
new task-oriented devices (such as quantum information 
transport, quantum cloning and simulation devices) con- 
stitute a new paradigm for quantum information process- 
ing (QIP), in which even a partial control over the system 
is enough to accomphsh the desired task. Of course, even 
for these devices a precise knowledge of their dynamics 
is an indispensable ingredient, thus spurring new inves- 
tigations and the development of new tools for coherent 
control. 

Here we study a physical system that presents a sim- 
pler quantum dynamics because of its one-dimensional 
geometry, while still retaining the full complexity of 
many-body dynamics. The system is a single crystal of 
fluorapatite (FAp) , where we focus on the spin degrees of 
freedom of the fluorine- 19 nuclei. This spin system has 
already been studied in the nuclear magnetic resonance 
(NMR) literature and proposed for quantum computa- 
tion. 

The goals of this study are to investigate the extent 
of our control on the system, when the control Hamil- 
tonians can only act collectively on all the spins; and 
to devise methods for extracting information on the sys- 
tem, when state tomography is not only inefficient, but 
also impossible because we can only detect the collec- 
tive magnetization. We do not consider here issues of 
decoherence mechanisms and rates, relying on the long 
coherence times of nuclear spins for performing all ex- 
periments in the coherent regime. Decoherence proper- 
ties of quantum many-body systems is however a very 
important area of study [U, @, S] j ^^'^ preliminary results 
regarding the system under investigation have already 
been presented 

In section II we first show how we can push the limits 
of our control when preparing a particular initial state, 
which is of interest for some tasks such as quantum in- 
formation transfer and quantum simulations [5| . Starting 
from the thermal equilibrium state, wc will show how to 



create a state where only the two spins at the ends of 
the chain are polarized, while all the other spins are in 
the identity state. Also, this state is the first step toward 
showing universal control on the system, as it has been 
proved that collective control on a qubit chain plus con- 
trol of the spins at the chain ends is enough for universal 
control of the system [1, 0] ■ In section III we present the 
characteristics of the physical system and the experimen- 
tal apparatus as well as the first indications that we suc- 
cessfully prepared the desired state. We then explain in 
section IV how we can extract enough information from 
the detection of the collective magnetization, to verify the 
creation of this desired initial state. The measurement is 
preceded by evolution under a propagator engineered to 
let emerge the many-body properties of the system. We 
use the technique of multiple quantum coherence (MQC) 
Q as a way to perform a partial tomography of the sys- 
tem and gain further information on its state. In particu- 
lar, we take advantage of the lower dimensionality of the 
system, which makes possible an analytical description 
of the evolution. We finally present experimental results 
in good agreement with the theoretical predictions, thus 
proving the preparation of the desired state. 



II. SELECTING THE SPINS AT THE 
EXTREMITIES OF THE LINEAR CHAIN 



There has been recently a remarkable interest in linear 
chains of spins [l^, [13, [13, [isll. as well as their simu- 
lation through atomic lattices |14| . The physical system 
we are interested in (a crystal of fluorapatite) can be 
modeled by a linear chain of spins-1/2 particles of a size 
(A^ > 10) such that many-body properties emerge. The 
system is put at room temperature in a large magnetic 
field and is addressed through radio-frequency (rf ) pulses. 
The internal Hamiltonian of the system in a frame rotat- 
ing at the larmor frequency is given by the secular part 
of the dipolar Hamiltonian: 



(1) 



while external control is provided by the rf Hamiltonian: 
Hrf ^ cj,f{t) [al cos (j){t) + al sin 0(i)] (2) 

i 

The control is therefore collective, that is, we can only 
address all the spins in the chain together: This limita- 
tion precludes universal control on the system. A recent 
scheme for QIP 0] showed however that adding to the 
collective control the ability to manipulate the spins at 
the chain ends enables universal quantum computation. 
A step toward achieving universal control would be to 
demonstrate at least a partial addressability of an indi- 
vidual spin. To prove the extent of our control on the 
system, we set up to prepare a particular quantum state, 
in which only the spins at the extremities of the chain are 
polarized: pg = 1 -|- e{a], + a^) (notice that this highly 
mixed state emerges from an ensemble average over many 
chains that are here considered as independent and equiv- 
alent systems). This state not only proves that we can 
break the internal symmetry of the system, but is also a 
useful state, for example for quantum information trans- 
port [H]. 

The spins at the extremities of the chain only have 
one nearest neighbor to which they are strongly coupled. 
This implies a different dynamics under the internal dipo- 
lar Hamiltonian with respect to the other spins. Taking 
advantage of this fact and complementing unitary con- 
trol with incoherent control will allow us to reach 
the desired state with a high enough fidelity. 

The spin system is initially at equilibrium, in the 
thermal state; in the high field, high temperature limit 
this can be very well approximated by the state peq = 

Ea^i "'z (notice that we ignore the large component pro- 
portional to the identity since it docs not evolve nor con- 
tribute to the signal). By rotating the thermal equilib- 
rium state to the transverse plane p{t = 0+) = J2k=i '^x^ 
the state is no longer an eigenstate of the internal Hamil- 
tonian and will therefore evolve. This evolution is usu- 
ally recorded in the NMR free induction decay (FID). It 
is known that the apparent decay of the magnetization 
masks a complex dynamics |17| . where non-detectable 



many-body states are created from the single-body mag- 
netization terms by the dipolar Hamiltonian couplings. 
Due to the fewer number of couplings to other spins, the 
first and last spins have a much slower dynamics (appar- 
ent decay) at short time. It is thus possible to select a 
particular time ti at which while the state of these two 
spins is still mainly <7x , sll the other spins have evolved to 
more complex multi-body states. From the commutator 
expansions of the unitary evolution, we can calculate the 
approximate coefficients of the polarization (cr^) terms 
for each spin as a function of time and therefore select 
the time at which a^^ Ri 0, V k ^ 1, N (because we are 
interested in the short time regime, an expansion to the 
gtft, Qj-jgj- gives an excellent approximation). Except for 
very short chains (3-4 spins), the polarization of all the 
spins fc > 2 is almost equal after this very short evolution 



and the optimal time is nearly independent of the num- 
ber of spins in the chain, therefore allowing us to choose 
the time ti even without knowing the precise (average) 
number of spins in a chain. 

A second tt/2 pulse will bring the magnetization of 
spins 1 and N to the longitudinal {z-)axis, so that the 
density matrix describing the system can be written as 
p{ti) = a{(j\ + a^) + p'. To select only the first two 
terms, which are the desired state, we can recur to a 
phase cycling scheme that selects only terms that com- 
mute with the total magnetization along z (X^^Li '^z) 
such as population terms. Unfortunately, we have not 
found a solution that also cancels out the zero-quantum 
terms (that is, components of the density matrix with 
total magnetic quantum number=0). Even with this lim- 
itation, the fidelity with the desired state is about 70%; 
the larger errors arc given by residual polarization on 
spins 2 and N-1 as well as correlated states of the form 



III. EXPERIMENTAL SETUP 

We performed experiments on a 300MHz Bruker 
Avance Spectrometer, with a home built probe tuned to 
282.4MHz for the observation of fluorine spins. The sam- 
ple was a single crystal of fluorapatite (Ca5(P04)3F). Ap- 
atites, either hydroxyapatites or fluorine containing ap- 
atites have been studied in NMR experiments be- 
cause of their particular geometry [23, Hlj and have also 
been proposed as a system to implement QIP [l^. The 
fluorapatite crystallizes in the hexagonal-dipyramidal 
crystal system, with cell dimensions a = 9.367A and 
c=6.884A, and two formula units per cell. The fluorine 
spins are arranged on linear chains along the c direction, 
with distance between two atoms d = 3.442A, and 6 ad- 
jacent parallel chains at the distance D ~ 9.367A. Since 
the dipolar couplings decrease with the cube of the dis- 
tance between spins, the spin system can be considered 
a quasi-lD system. In particular, by orienting the c axis 
of the crystal along the z direction, the ratio of the cross- 
chain to in-chain dipolar coupling is: 



6x 1 <P 

= 0.0248. 

kn 2 



(3) 



Although this crystal has a very good purity, as testi- 
fied by long decoherence times (Ti ~ 200s), the chains 
are interrupted by defects and vacancies, that limit their 
length. If we consider the short-time evolution only, the 
spins of different chains do not interact and behave as in- 
dependent systems, while if we let the system evolve for 
a time long compared to the cross-chain interactions, the 
approximation of independent systems will fail. We will 
include the effects of other chains and the distribution 
in chain length into the environmental decoherence and 
only consider a single chain as the system of interest . 

We applied the pulse sequence presented in the prcvi- 



2 



ous section: 



IV. MULTIPLE QUANTUM COHERENCE AS A 
STATE MEASUREMENT TOOL. 



where the pulse axis a was phase cycled through y and 
X to cancel the non zero-quantum terms. The optimal 
time tl for the dipolar coupling strength of fluorapatitc 
spin chains is given by ti = 30.3/is. 

In order to assess qualitatively the results of this se- 
quence, we measured a spectrum of the system, after the 
preparation of the desired initial state and compared it 
to the thermal equilibrium spectrum. In particular, it is 
important to observe the state at short time, when the 
effects of the desired initial state should be stronger; the 
solid echo technique [1^ was thus used, to avoid the dead 
time imposed by the electronics and the pulse ring-down 
[2^ . While this kind of measurement does not give a 
definite answer to the question whether the polarization 
has been concentrated on the extremities only, the qual- 
itative differences in the spectra measured are encourag- 
ing. In particular, we observe that resolution of the three 
peaks observed in the thermal state is much better (we 
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FIG. 1: Comparison of the spectra when the polarization 
is retained by all spins in the chain (dashed line) and for the 
excitation of the extremities only (solid line). The FWHM 
is « 19kHz and the distance between peaks ~ 8kHz. The 
experimental data were obtained with a modified version of 
sequence Q and a Solid Echo read-out. (The sequence was 
equivalent to the one presented, but an extra pulse gave exper- 
imentally a cleaner spectrum when combined to the read-out 
sequence). The wait time ti was ti = .5/xs for the all spin 
spectrum and 30.3^is for the chain ends excitation. 

obtain narrower lines) as expected from a state in which 
fewer couplings are available, and while the spectrum of 
the thermal state presents the characteristic of a triplet 
[Tot as expected for one spin coupling strongly to 
two equivalent spins (as is the case if we consider nearest 
neighbors only) the spectrum for our state is closer to 
a doublet, reflecting the fact that the spins at the chain 
extremities interact strongly with only one spin. 



Since the detection of magnetic resonance restricts the 
obscrvables to the collective transverse magnetization, it 
is not possible to reconstruct the state of the system by 
quantum state tomography [2^ . To assess the efficiency 
of the state preparation scheme, we thus need a read- 
out scheme that reveals the signature of the particular 
initial state prepared. By studying a more complex dy- 
namics than the free induction decay, we can obtain more 
accurate insight into the state created. In particular, a 
sensitive probe of the dynamics of a correlated many-spin 
system is the creation and evolution of quantum coher- 
ences. 



A. Fermion operator solution to MQC dynamics 

Quantum coherence refers to a state of a physical sys- 
tem where the phase differences among the various con- 
stituents of the wave function can lead to interferences. 
In particular, quantum coherences often refer to a many- 
body system, whose parties interact and therefore show 
a correlation, a well defined phase relationship. In NMR, 
coherences between two or more spins are usually called 
multiple quantum coherences. When the system is quan- 
tized along the z axis, so that the Zeeman magnetic mo- 
ment along z is a good quantum number, a quantum co- 
herence of order n is defined as the transition between two 
states \mi) and |to2), such that the difference of the mag- 
netic moment along z of these states is n: nii — TO2 oc n. 
Multiple quantum coherences of order n usually describe 
states like \m2){mi\, or elements in the density matrix 
that correspond to a transition between these two states 
[HjII^. The state |m2)(mi| is also called a coherence of 
order n. Diagonal density matrix terms are called popu- 
lations and do not properly describe a coherence. 

The observation of multiple quantum coherences in 
NMR started in the mid 1970s, as a method for unrav- 
eling complex spectra, by filtering transitions based on 
the coherence order involved [23|. More recently, these 
states have been studied with respect to their decay time 
and effects of decoherence on them [l|, 0] , because of their 
connection to entangled states. Multiple quantum coher- 
ences are also the basis of proposed schemes for the full 
characterization of complex many-spin states (28j : iden- 
tifying the different coherence orders occurring in a state 
is the first step toward quantum state tomography. 

Coherences can be created by the interplay of rf pulses 
and free evolution periods, during which interactions 
among spins occur. Multiple pulses sequences create 
Hamiltonians that can raise the coherence order of the 
system. One widely used sequence [1, [2^ , for example. 
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creates the double-quantum Hamiltonian 



^y„y 



(5) 

which raises the coherence orders in step of two (thus, 
starting from the thermal state, only even coherence or- 
der will be created) . The growth of coherence orders has 
been the matter of many investigations [13, [lOl , because 
of its relationship to the system's geometrical structure. 
Because of its complexity, only stochastic models, based 
on the probability of occupation of different coherent 
states, and semiclassical models like the hopping model 
are available to describe this dynamics. If we restrict 
the evolution to one-dimensional systems, the growth of 
coherence is slowed down by the fewer couplings among 
spins; if we further assume that only nearest-neighbor 
couplings are present, the evolution under the double- 
quantum Hamiltonian turns out to be exactly solvable 
jsil . [3^ . Experimentally, the nearest-neighbor (nn) ap- 
proximation is accurate for short times, while for longer 
times, weaker couplings start to produce appreciable cor- 
rections. 

The most important characteristic of 1-D MQC ex- 
periments are the oscillations between zero- and double- 
quantum coherences at short times. It is this restric- 
tion of the accessible Hilbert space (that is exact at any 
time in the nn approximation) that makes the problem 
analytically tractable. As we will show in this section, 
these oscillations turn out to be a signature of the initial 
state, so that they can help us confirm experimentally 
the preparation of the initial state desired. For compari- 
son, we derive analytically (and then measure experimen- 
tally) not only the MQC intensities for the desired state 
/9e(0) = <tI + , but also for the thermal equilibrium 

state p{0) = J2j=i'^i- 

The analytical result is obtained by mapping the spin 
system to spinless fermion operators (see the appendix): 



c,= n 

k=l 



(6) 



The diagonalization of the double-quantum Hamiltonian 
[31I [3^.l33t is accomplished by a Bogoliubov canonical 
transformation to the operators dk [3J]: 



1 ^ 

p^^sin(Kj) (7fcrffc+rfLfc), 



n = (7 

N + l ^ ' 



where 7^ = sgn(/c). The DQ-Hamiltonian is then diago- 
nalized to: 



Udq = -26 ^ cos At {d\dk + d^fcrf-fc - 1) (8) 



fc=i 



We now express the two initial states in terms of Bo- 



goliubov operators d^. The thermal state has a particu- 
larly compact expression. Using first the spin to fermion 
mapping, we have 



N N , 



TV 



N 



= -5](d,d_,-4dl,) 



fe=i 



(9a) 



(9b) 



where in writing (j9b[) we used orthogonality relationships 
for trigonometric functions (see Eq. ([3T|) and (|32p in the 
appendix) . 

Consider now the initial state Pe(0) = + cr|^). 

Because we are not summing over all spins, it is no longer 
possible to use the orthogonality relationships as in (j9b[) . 
This results in more cumbersome double sums: 

Pe(0) = \{al cr^) = 1 - (cjci + cjvCAr) (10a) 

= 1 - E ^^'^(7^4 + d^uMnudu + 4j, (10b) 

k,h 

where to simplify the notation we set S-^h = 
sin(K) sin(77) -|- sin {N sin (Nrj). 

The system evolves under the double quantum Hamil- 
tonian with a dynamics described by the propagator 
U{t) = exp {^ibj^k cosK{dldk + d^_f.d-k))- We now de- 
fine the eigenphases V'fe = 2btcoaK. The thermal state 
evolution is easily calculated to be: 



Pit) 



1 



dkd-kc 



-2'iV'fc 



(11) 



In order to separate contributions from different co- 
herence orders, we transform back to fermion opera- 
tors. To simplify the calculations, we use the operators 
a-k — (7fedfe + (iLfc)/'\/2 that represent the same coherence 
order as the c^'s operators (the dk operators correspond 
instead to different coherence orders, since they are com- 
binations of lowering and raising operators) . With some 
algebraic manipulations, we have 



~ ^ E 7fc sin IV-fel (alalfc + afcO-fc) 



(12) 



The intensities J„ of each 



quantum coher- 



ence as measured in MQC experiment is given by 
Tr [((7z3qp(0)C/|,q)(")(C/z,q E^^.'^dq)^"^] • Since in this 
case p{0) = Y.a^,we have J„ = Tr . We eval- 

uate the trace of the fermion operators Uk , aj. in their cor- 
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responding occupational number representation, so that 
only terms like a^Uk are diagonal and contribute to the 
trace (see Eq. psp and in the appendix). The nor- 
malized MQC intensities for zero and double quantum 
arc finally given by [s^ : 

Jo = ^ (^^^ ^'^^ (13a) 

k 

J2 = sin^ (Abt cos k) 




(13b) 



MQC excitation time (ms) 



The evolution of the polarization on the extremities of 
the chain is : 



Again, we transform to fermion operators to distin- 
guish contributions from different coherence orders: 



Peit) 



i 1 



2N- 



k.h 



-1 - 



N 



~rj X! ^kh[alah cos {ipk + iph) + 5k,h sin V'fe^] 

k.h 



(15) 

As expected, also starting from the non collective ini- 
tial state only zero and double quantum coherences are 
developed (taking into account nearest neighbor cou- 
plings only). 

The zero quantum contribution to the signal is Jg = 



Tr 



(with p(°)(i) from eq. (HH)): 



alak - \ 



Jq = -EfcCOSi/jfeTr 
+ (IvTT) X! '^"^ ^' ^'^^ + V'/jTr 



k,hj 

hpT) X! cos -0; sin Tr 



a\aha\ai 



k,l 



t -L 



(16) 

Using the trace and trigonometric identities in the ap- 
pendix, we find the normalized zero-quantum intensity: 

Jo = ^i^' ('^) cos At) (17) 



The double quantum intensity J| is given by 
JTm)T,kM,iSkhSini;i 



Jl =Tr 



X i sin [tpk + V'/i)Tr 



(18) 



FIG. 2: Zero- and double-quantum intensities as a func- 
tion of the evolution time under the double-quantum Hamil- 
tonian. Nearest-neighbor couplings only are assumed, with 
equal strength as given by the fitting to experimental data 
(see Fig. (O). In particular notice the clear difi'erences in 
the behavior for the two initial states. Also the even-odd 
spin number dependence of the MQC intensities is interest- 
ing: while this tends to go to zero for large number of spins 
in the collective initial state case, this difference is observed 
even for very large number of spins for the other initial state. 



Y '^fc'' ^^^"^ (^'^ + V^'0(ala-fc - a-kan) which yields: 



J2 = W+T) Sfc si'^^ (^) '^"^^ (4^* cos k) 



(19) 



These more complex expressions lead to a very differ- 
ent behavior of the coherence intensities as shown in Fig. 
(I2|), so that it is possible to distinguish even experimen- 
tally what was the initial condition of the system. 



B. Experimental results 

We applied the pulse sequence (|4]) followed by a MQC- 
experiment sequence. In particular, we used the 16-pulse 
sequence p9j (except for the 3 shorter time- values, where 
the 8-pulse sequence was used Q), and a phase cycling 
with increments of 1^9 = 27r/8 to select up to the 4*'' 
coherence order, by repeating the experiment 16 times. 

The time delay between pulses in the MQC sequence 
was varied from 2iis to 6.5/xs, to increase the excitation 
time, as well as the number of repetitions of the sequence 
itself (1 or 2 loops), so that the evolution of the quantum 
coherences were studied between the times of 37.6/is to 
354.4/is. We compare the results obtained with the evo- 
lution for the thermal equilibrium as initial state. In or- 
der to take into account the effects of imperfections in the 
pulse sequence, we applied the pulse sequence as in Q 
also to obtain the thermal state, with a very short ti time 
{ti = O.bfis). In figure ([3|) we show the dynamics of the 
zero and double quantum intensities, normalized to have 
sum = 1 to take into account the signal decay for longer 
excitation times. We notice that the four-quantum coher- 
ence intensity is as low as the baseline, indicating that 
the time scale is short enough for the nearest-neighbor 
approximation to be valid to a good extent (remember 
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.05 1 .15 .2 .25 .3 .35 .05 1 .15 .2 .25 .3 .35 

MQC excitation time (ms) MQC excitation time (ms) 



FIG. 3: Experimental results. Left: the initial state is the collective thermal state X^^cr^. The experimental points have 
been fitted (dashed line) to the theoretical curves for nearest neighbor coupling only, with the dipolar coupling as fitting 
parameter. The number of spins was varied to find the best fit, which results to be A'" = 11 spins. Right: MQC intensities for 
the initial state po = o"i +ct^ . Also plotted are the theory predictions for the same dipolar coupling and 11 spins (solid line) or 
10 spins (dashed line). The experimental behavior observed (a constant behavior also for longer time) indicates that a simple 
model with 11 spin chains is not an accurate description at longer times (see text). 



that the nn-approximation predicts that only zero- and 
double-quantum coherences are excited). 

The experimental results for the MQC oscillations 
starting from the thermal state have been fitted to the 
theoretical curve (|13p for a single chain with nearest 
neighbor couplings only. The dipolar strength and the 
number of spin were the fitting parameters. The results 
of this fitting were used to plot the theoretical curve for 
the chain ends initial state [Eqn. pT]) and (HH)]. The 
concordance of the theoretical predictions with the ex- 
perimental data is very good, even if the state created 
contains residual zero-quantum terms. 

The best fitting of the experimental data in Fig. 3 to 
was found for a chain length of 11 spins. This re- 
sult must be taken with caution, since the chain length 
influences only the long time behavior of the MQC inten- 
sities, where other factors not taken into account in the 
simple analytical model start to play an important role. 
In particular, the experimental results for the chain ends 
indicate that there are either longer chains or a distribu- 
tion of short chains, with odd and even number of spins. 
Although the reason for the best fit at 11 spins could 
be simply due to a distribution of spin chains around a 
mean of 11 spins, the low impurities content of the crystal 
studied is not in agreement with this finding. Wc believe 
that a more plausible reason is the break down of the 
model of isolated nearest-neighbor coupled spin chains. 
Even if next-nearest neighbor couplings and couplings to 
adjacent chains are not strong enough yet to create four- 
quantum coherences, they still modify the intensities of 
the zero- and double-quantum coherences. Further ex- 
periments would be needed to distinguish between these 
hypotheses. 



V. CONCLUSIONS 

We have studied a naturally occurring physical sys- 
tem, a single crystal of FAp that presents linear chains of 
spins- 1/2 particles. Since the physical characteristics of 
the system and the experimental apparatus do not pro- 
vide universal control on the quantum spin system, we 
propose to use this system not as a candidate quantum 
computer, but as a specific task-oriented QIP device, for 
example for quantum state transport or simulations. In 
particular, we have devised a scheme, combining unitary 
and non-unitary control, for creating a particular state 
that breaks the natural symmetry of the system. This 
state will allow us to study properties of state transport 
along the chain as well as being an interesting initial 
state for simulations. Furthermore, the preparation of 
this state is the first step toward universal control on the 
system, since full control on the spins at the chain ends 
(in addition to the collective control over all other spins) 
ensures universality. 

In addition, we have investigated a tool for acquiring 
a deeper knowledge of the state and dynamics of the sys- 
tem, given the limitations in the read-out procedures. 
Multiple quantum coherences allow us to gather more in- 
formation on the multi-body aspects of the system than 
simple direct observation of the collective polarization. 
In particular, we used analytical solutions in the limit of 
nearest-neighbor couplings to interpret the experimental 
results which confirm the preparation of the desired state. 

In conclusion, we have shown how even a quantum 
system without universal control can be used to study 
physical problems of interest in condensed matter theory 
and quantum information science. 
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VI. APPENDIX: FERMION OPERATORS 

Spin operators of the Pauli group can be mapped to 
fermion operators, obeying the well-known anticommuta- 
tion relationships. This mapping is useful in describing 
the dynamics of various ID models, since some Hamilto- 
nians can then be diagonalized analytically. In the fol- 
lowing we describe a particular mapping that is suited 
for describing the creation of MQC. 

A mapping from spin to fermion operators goes back 
to Jordan and Wigner [3^, who first transformed quan- 
tum spin S = 1/2 operators, which commute at differ- 
ent lattice sites, into operators obeying a Clifford alge- 
bra (fermions). This transformation was used to map 
the one-dimensional Ising model into a spinless fermion 
model, which is exactly solvable. The Jordan- Wigner 
transformation has been recently generalized to the cases 
of arbitrary spin S [s^, [13] and to 2D spin systems [1^ . 

Given a set of spin-i operators cr", each defined at a 
lattice site j, they obey the commutation relationship: 



cr., ,0- 



Sj,kio-' 



(20) 



where {a,/?, 7} = {x,y,z} and cyclic permutations of 
these indexes. The raising and lowering operators cr^ = 
((Tj±icrJ)/2 obey mixed commutation and anticommuta- 
tion relationships, which are not preserved under a uni- 
tary transformation. To diagonalize the spin Hamilto- 
nian we must thus use the Jordan- Wigner transforma- 
tions to map these operators to fermion operators c,,- , , 
obeying the canonical anticommutation relationships: 



{c],Ck} - S,,k, {c„ck} - {c], 4} = 0, (21) 

where we adopted the notation { , } for anticommu- 
tators. A basis for the Hilbcrt space of these opera- 
tors is given by the occupation number representation 
\n) = 712, JT-Af), where Uj ~ {0,1} is the occupa- 
tion number at site j. The state \n) can be obtained from 
the vacuum state by: 



11(4)"^ 



in) = I I (c^)-'\vac) 



(22) 



Then, the action of the fermion operators on such states 
is given by: 

c,N = (°V , ^'"^■=° (23) 

' [-{-l) i\n'), otherwise ^ ' 

where \n') is the vector resulting when the j*'* entry of 



is changed to and s" = X]i=i Analogously: 
0, 



cl\n) 



if ?ij = 1 



-(— l)^j |n'), otherwise 



(24) 



where \n') is the vector resulting when the j*^ entry of 
\n) is changed to 1. 

The mapping from spin to fermion operators can be ex- 
pressed in several ways, the most intuitive being based on 
identifying every basis vector |n) in the occupation num- 
ber representation basis to the corresponding \n) basis 
vector in the computational basis for the spin operator 
Hilbcrt space. Imposing this one to one correspondence 
on basis states and taking into account the respective ac- 
tions of spin and fermion operators on their basis vectors, 
one obtains the mapping: 



k=l 

J _ _ TT ^+ 



n 



k=l 

Notice also that erf 



'=1 (25) 

n(i-24c.)ct 

k=l 



l-2c]c,. 



Consider now the double quantum Hamiltonian with 
equal couplings restricted to the nearest neighbor spins : 



TV 



3 = 1 



(26) 



We can express it in terms of the fermion operators as: 

-Hdq = & En(l - 24q)c]+i - 24q)c] + h.c. 
j=i 1=1 1=1 

N 

(27) 

where we introduced the vector — [c^c] and the matrix 
B: 



B 



B 

-B 



B 



(28) 



Notice that even if the mapping to fermion operators 
is non-local, this quadratic Hamiltonian is mapped to 
a local Hamiltonian, when one considers only nearest 
neighbor couplings. The matrix B is orthogonal, and 
can therefore be put into diagonal form, with eigenvalues 

cos (^^^^ , k = 1, 2, ... TV and eigenvectors = [d^d]. 

The diagonalization is a canonical unitary (orthogonal) 
transformation (a Bogoliubov transformation) to a di- 
agonal basis satisfying the same anti-commutation rela- 
tionships [33| . To simplify the diagonalization we need 
to introduce also the negative modes operators dk and 
correspondingly the fermion operators Cj (with j < 0), 
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by extending the sum in (|27p from —N to N. 

We can express the fermion operators cj in terms of 
the bogoHubov operators dk by the following orthogonal 
relationship: 



1 ^ 

==^sin(Kj)(7fc4+dLfc)> 1^ 



irk 
N+l 



(29) 



where "fk = sign(fc), ensuring the transformation is or- 
thogonal (this transformation is obtained by imposing 
that the resulting DQ-Hamiltonian is diagonal). Substi- 
tuting expression (|^^ in equation (P7|) wc have: 



N 



k.h=l 
1 ^ 

^TVTT ^ sm{Kj)sin[j]{j + 1)] 

j = -N 

(30) 

(with Tj = nh/{N + 1)). By using the following orthogo- 
nality relationships: 



N 



1 ^ 

— ^ sin(fcj) sin(/ij) = {5k^h - 4 -h) (31) 



j = -N 



and: 



AT 



sin(fcj) cos(/ij) = 



(32) 



we can simplify the Hamiltonian pO|) : 



T^CQ = -25^ COS At {dldk + rfLfcrf-fc - 1) (33) 
fc=i 

We note that under this diagonal Hamiltonian, the 
operators dk evolve as: dk{i) — e^*'^°<'*(ifce'^°Q* = 

^ — 2ib cos 

o Uk • 



For convenience, we introduce also fermion operators 
Ofc that are useful in expressing the density matrix: 



ak = {ikdk + dl_^)/V2 



(34) 



The following identities for the trace of fermion opera- 
tors have been used in the text to calculate the zero- and 
double-quantum intensities : 



Tr 



Tr 



alak 



a\aha[,aki 



Tr 



a\a^_f,a^k'ah' 



2^-2 forfc^fc'; (35) 
2^-1 for k^k'. 



2^-2 for fc = fc',/i = /i'; 
-2^-2 foj, ^ _.h\h = -k'; 
otherwise 

(36) 
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